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ABSTRACT 

Using time-dependent linear perturbation theory, we evaluate the dynamical 
friction force on a massive perturber M p traveling at velocity V through a 
uniform gaseous medium of density p and sound speed c s . This drag force 
acts in the direction —V, and arises from the gravitational attraction between 
the perturber and its wake in the ambient medium. For supersonic motion 
(M. = V/c s > 1), the enhanced-density wake is confined to the Mach cone 
trailing the perturber; for subsonic motion (Ai < 1), the wake is confined to a 
sphere of radius c s t centered a distance Vt behind the perturber. Inside the wake, 
surfaces of constant density are hyperboloids or oblate spheroids for supersonic 
or subsonic perturbers, respectively, with the density maximal nearest the 
perturber. The dynamical drag force has the form Fjf = —I x 4n(GM p ) 2 po/V 2 . 
We evaluate / analytically; its limits are / — > Ai 3 /3 for M. « 1 and 
/ — > \n(Vt/r min ) for M. » 1. We compare our results to the Chandrasekhar 
formula for dynamical friction in a collisionless medium, noting that the gaseous 
drag is generally more efficient when Ai > 1 but less efficient when A4 < 1. To 
allow simple estimates of orbit evolution in a gaseous protogalaxy or proto-star 
cluster, we use our formulae to evaluate the decay times of a (supersonic) 
perturber on a near-circular orbit in an isothermal p oc r~ 2 halo, and of a 
(subsonic) perturber on a near-circular orbit in a constant-density core. We 
also mention the relevance of our calculations to protoplanet migration in a 
circumstellar nebula. 



1. Introduction 

The process of dynamical friction, defined as momentum loss by a massive moving 
object due to its gravitational interaction with its own gravitationally-induced wake, arises 
in many astronomical systems. Examples of systems in which such effects are well known 
to be important include stars in clusters or galaxies, galaxies in galaxy clusters, and binary 
star cores in the common envelope phase of evolution. In the first two examples, the 
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surrounding background medium in general consists of a combination of collisionless matter 
(stars, galaxies, dark matter) and gas, while in the third example the surrounding medium 
is entirely gaseous. As a corollary to the dynamical friction process, in all these cases the 
background medium is heated at an equal and opposite rate to the energy lost by the 
perturber. 

The analytic theory for the gravitational drag in collisionless systems was developed 
by phandrasekhar (1943)| , and over the decades since has enjoyed widespread theoretical 
application, extensive verification by numerical experiments, and well-documented 
embodiment in observed astronomical systems. The variety of important consequences of 
gravitational drag in collisionless astronomical systems includes mass segregation in star 
clusters, sinking satellites in dark matter galaxy halos, orbital decay of binary supermassive 
black holes after galaxy mergers, etc.; see e.g. [Binney fc Tremaine (1987)] . 

Less well-developed is the corresponding theory of dynamical friction in a gaseous 
(i.e., collisional) medium.§ For supersonic motion, analytic linear-theory estimates of the 
gravitational drag under assumption of a steady state were obtained by pokuchaev (1964) 



Ruderman fc Spiegel ( 197TJ] , and [Rephaeli fc Salpeter (1980]| . The resulting estimates for 



the drag force in the steady supersonic case take the form 

An(GM p ) 2 p 



F ss = " y. In 

where r max and r min correspond respectively to the effective linear sizes of the surrounding 
medium and the perturbing object, similar to the drag formula obtained for a collisionless 
medium. Although there is some ambiguity in the definition of r max and r min , the estimate 
(|l|) appears consistent with calculations of the gravitational drag that are obtained as a 
by-product of numerical hydrodynamic investigations focused on the Bondi-Hoyle-Lyttleton 
accretion problem - see, e.g. |Shima et al (1985)| , |Shankar et al (1993)| , and Ruffert and 



collaborators (see Rufiert (1996) and references therein). 



For steady-state, subsonic motion of the perturber, the front-back symmetry of the 
perturbed density distribution about the perturber led [Rephaeli fc Salpeter (19801 to 



argue that gravitational drag is absent in the subsonic, inviscid case. Considering that the 
drag force for the supersonic case increases proportional to V~ 2 with decreasing perturber 
speed, V, it seems counterintuitive for the dynamical drag to become exactly zero when 
V becomes infinitesimally smaller than the sound speed. In this paper, we reconsider the 



instead, most studies of the gravitational interaction between a moving massive body and the surrounding 
gaseous medium have focused on the problem of accretion, following on the analysis of Hoyle and Lyttlcton 



(1939) and Bondi and Hoyle (1944), and the early numerical work of Hunt (1971 
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linear-theory drag as a time-dependent rather than steady-state problem, and arrive instead 
at a nonzero value for the dynamical friction for the subsonic case, while still verifying that 
the drag force is maximized for perturbers with V ~ c s . In §2, we derive results for the 
perturbed density distributions created by, and dynamical drag on, a massive perturber 
on a constant-velocity trajectory through a uniform, infinite medium. In §3 we relate 
our results to the classical (collisionless) dynamical friction formula, and briefly consider 
applications to a few astronomical systems. 



2. Analysis 

2.1. Wave Equation and Formal Solution using Green Functions 

We begin with the linearized equations for the perturbed density p = p [l + a(x, t)) 
and velocity v = c s /3(x, t) of an adiabatic gaseous medium that is subject to an external 
gravitational potential $ ex t(x, t): 

^ + ^ = ° (2) 
c s at 



and 

c* dt c 



+ V« = -^V$ext, (3) 



where po is the unperturbed density, c s = (dPo/dpo) 1 ^ 2 is the sound speed, and the 
perturbation amplitudes a, \/3\ « 1. By substituting equation (0) in the divergence of 
equation (|3|), we have 

V 2 « - ~ = -^V 2 $ e „ = -47r/(x, t) (4) 
where p ex4 (x, t) = c^/(x, t) is the mass density of the perturber. 

In the absence of any disturbance prior to the action of $ ex t, the solution to equation 
(|) is found using the retarded Green function G + for the 3D wave equation QJackson (1975) 
§6.6) as 

a(x, t) = ( I M /V-{t-\*-Ah,)\f(*,<) (5) 



2.2. Perturbed Density Distributions for Constant- Velocity Perturbers 



We now specialize to the case of a point mass M p on a straight-line trajectory with 
velocity Vz, passing at time t = through x = 0. If T~C{t) describes the time over which the 
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perturber is active, then /(x, £) = (GM p /c 2 )5(z — Vt)6(x)6(y)H(t) . We define s = z — Vt as 
the distance along the line of motion relative to the perturber, w = z' — z, and M. = V/c s as 
the Mach number; w > 0(< 0) corresponds to the perturbation from a backward (forward) 
- propagating wave. Then 

„U A - GM p f°° j J-t + s + M(R 2 + w 2 y/ 2 }H((w + z)/V) 

alX ' tj " C 2 LJ W (^ + ^1/2 • W 

Here R = (x 2 + y 2 ) 1 ^ 2 is the cylindrical radius. 

To evaluate the integral (|6]), one expands the argument of the 5— function about its 
possible roots 

s±M[s 2 + R 2 (l-M 2 )W 2 

for M. < 1, only w + is a valid root, whereas for M. > 1 both roots are valid as long as s < 
and \s\/R > (M 2 — l) 1 / 2 , and neither otherwise. Using 8((w — w±)A) = 5(w — for 
A = 1 + Mw±/(R 2 + w 2 .) 1 / 2 and substituting the solution (0), the result is 

CM I r 2 

For a purely steady solution over all time, 7i is unity for all arguments, with the result 

that 

CM lc 2 f 1 \IM<1 

as = W + R 2 {l P _M2 )Y /2 >< 2 if M > 1 and s/R < -{M 2 - if' 2 (9) 

I otherwise 

This confirms the analytic, linear-theory results of previous authors for completely steady 
flow created by a point mass on a straight-line, constant-speed trajectory: (a) A subsonic 
perturber generates a density distribution centered on the perturber at s = 0, with contours 
of constant density corresponding to similar ellipses in the s — R plane with eccentricity 
e = Ai, and the short axis along the line of motion of the perturber (i.e. the 3D density 
distribution consists of concentric similar oblate spheroids). The density contrast p 1 /p = a 
is unity for an elliptical section with semiminor axis GM p /c 2 ; the linearization is sensible 
only outside of this elliptical section. This density distribution is a generalization of 
the far-field limit of the hydrostatic envelope p/po = exp(GM p /(c 2 r)) that surrounds a 
stationary perturber.! (b) A supersonic perturber generates a density wake only within the 



2 We note that for subsonic flow, the far-field density enhancement from equation (||) is within 10% of 
Hunt's (1971) M. = 0.6 numerical solution outside of r = GM p /c 2 s . Also, the rms anisotropy of 7% from 
the far-field steady solution (||) with M = 0.6 is comparable to the 5% mean anisotropy cited for Hunt's 
near-field numerical solution. 
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rear Mach cone of half-opening angle sin# = 1/Ai defined by s/R < — (Ai 2 — l) 1 ^ 2 ; the 
surfaces of constant density within the wake correspond to hyperbolae in the s — R plane, 
with eccentricity e = Ai. 

Now consider the case where the perturber is "turned on" at t = 0, so that TC is a 
Heaviside function. For a subsonic perturber Ai < 1 the only root is w + (eq. 0). Algebraic 
calculation shows that z + w + > when R 2 + z 2 < (c s t) 2 , for Ai < 1. Thus, the region 
of perturbation is the sphere centered on the original position of the perturber, within 
which a sound wave has traveled in time t. Within this region of perturbation, the density 
distribution has reached the value given by the steady solution (eq. ||); outside the causal 
region for sound waves, the density remains unperturbed. 

For a supersonic perturber, any density disturbance must be confined within the rear 
Mach cone s/R < —(Ai 2 — I) 1 / 2 . Algebraic calculation shows that within the sphere 
R 2 + z 2 < (c s t) 2 , z + w + > and z + W- < 0; hence only w+ contributes in equation (§). 
Within the Mach cone and to the right (R < \z — Mc s t\/(M 2 - 1) 1/2 , z > c s t/M) of this 
sphere (R 2 + z 2 > (c s t) 2 ), both z + w + and z + w_ are real and positive, hence contribute 
in eq. (|8]). 

The results for the perturbed density for this finite-time perturbation are summarized 



as 



follows: i 



, , 1 if R 2 + z 2 < (c s t) 2 

a(t)= D2 . p/ him72 x 2 XM>l,R 2 + z 2 >(c s t) 2 , s/ R < -(M 2 - if' 2 , and z>c s t/M 
[s z + R z {± - M z ) Y IZ J _. 

I otherwise 

(10) 

The region of perturbed density has the shape of a loaded ice-cream-cone dragged by its 
point by the perturber M p ; only w + contributes in the ice-cream region ("region 1"), while 
both w± contribute in the cone region ( "region 2" ) . The cone shrinks in size as Ai decreases 
for Ai > 1, and is nonexistant for Ai < 1. Figures 1 and 2 show examples of the perturbed 
density distributions for subsonic and supersonic perturbers, respectively. Because of the 
linear-theory assumptions made at the outset, we note that equation (|10"D is properly valid 
only for a « 1, i.e. s 2 + R 2 (l - M 2 ) » (GM p /c 2 s ) 2 . 



3 These results appear to have been obtained previously by Just fc Kegel (1990) via an alternative 
mathematical formalism; the two calculations serve as independent checks of the formulae. 
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Fig. 1. — Density perturbation profiles for subsonic perturbers with M. = 0.1,0.5,0.75,0.9 
(as indicated in upper right). Contours show isosurfaces of log (a) = log(a) — log[GM p / (tc^)], 
in intervals of 0.1. Density increases toward the perturber; the heavy contour indicates the 
surface with a — 1. The + symbol indicates the initial position of the perturber. 




Fig. 2. — Density perturbation profiles for supersonic perturbers with M = 1.01, 1.5,2, 10 
(as indicated in upper right). Contours show isosurfaces of log (a) = log (a) — log[GM p / (tc^)], 
in intervals of 0.1. Density increases toward the perturber at the apex of the Mach cone. 
The heavy contour indicates the surface with a — 1. There is a density jump with 
Alog(a) = log(2) = 0.301 at the surface R 2 + z 2 = (c s t) 2 . The + symbol indicates the 
initial position of the perturber. 
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2.3. Gravitational Drag Formulae 

To compute the gravitational drag on the perturber, we need to evaluate the 
gravitational force between the perturber and its wake, 

F df = 2nGM pPo J J dsdRR {g2 a ^ 2)3/2 . (11) 

To perform the volume integral, it is convenient to transform to a spherical polar coordinate 
system (r, 9) centered on the massive perturber, with R = rsin# and s = rcos8. Defining 
ji = cos 9, x = r/c s t, we have 

F(r = _„, * s ±«* (12) 



V 2 



1 r dx r jiM 2 Su 



where 

I = ~2~J ~J dfX (l-M 2 + fitM 2 ) 1 ' 2 ' (13) 
Here Su represents the sum in equation (§). 

For the purely steady-state density perturbation given in equation (||), in the subsonic 
case S-h = 1 everywhere in space and front-back antisymmetry in the angle integral (due 
to the symmetric spheroidal density distribution) argues that there is zero net force on 
the perturber. For the steady-state, supersonic case, Su = 2 for all \i between —1 and 
Hm = —\/M. 2 — 1/M. (the boundary of the Mach cone), so that I — J dx/x = hi(r ma2 ./r m j n ). 
This is consistent with previous results, and yields an identical formula to that representing 
the dynamical friction force in a collisionless medium when V is much larger than the 
background particle velocity dispersion (e.g. [Binney fc Tremaine [1957] , eq. 7-18). 



For the finite-time case, based on the perturbed density distribution in equation fllPl), 
there is nonzero contribution to the integral I only from a finite region in space: region 
1, in which S-j-i = 1, and region 2, in which Sn = 2. Region 2 is nonexistant for subsonic 
perturbers. For the subsonic case, the perturber is surrounded by a concentric distribution 
of similar ellipsoids, and is displaced by Vt forward from the center of the sonic sphere 
(radius c s t) surrounding its initial position. The nearby, complete ellipsoids exert no net 
force on the perturber, but the larger, cut-off ones with semiminor axes between (c s — V)t 
and (c s + V)t lag behind the perturber (see figure 1) and exert a gravitational drag. Thus, 
the radial integral in equation (0) has upper/lower limits x = 1 ± A4, and the angular 
integral has limits /i = — 1,/ic with ft c = (1 — M 2 — x 2 )/{2xM). The integrals are 
straightforward; the result is 

I subsonic =\ ln ( tZTJ^i ) ~ M - ( 14 ) 
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The implicit assumptions in deriving this equation are that (c s — V)t exceeds the effective 
size of the perturber (r min ), and that (c s + V)t is smaller than the effective size of the 
surrounding gaseous medium (r max ). Under these conditions, the dynamical drag is 
time-independent, and non-zero. The steady-state result that zero net force results from 
the front-back symmetry of the density distribution is misleading; because of the long-range 
nature of the Coulomb potential, the total drag force at any time depends on the unchanging 
ratio (c s t + vt)/(c s t — vt) = (1 + M.)/ (1 — M.) between the semiminor axes of the furthest 
and closest perturbing partial spheroids. The gravitational drag is always dominated by 
the far-field, and at any time the perturber is located ahead of center of the sonic sphere. 
Physically, we can associate the energy loss arising from the drag force with the rate at 
which the expanding sound wave does work on the background medium. In the limit of a 
very slow perturber M. « 1, I su bsonic —* M. 3 /3, so that the drag force is proportional to 
the perturber's speed V. 

For the supersonic case, the whole of the perturbed density distribution lags the 
perturber. The angular integration limits are fi = — 1, Hm f° r x — r min/(c s t) to V ' M? — 1, 
and \i = —1,(jlc for x = V ' M 2 — 1 to M. + 1; Sn takes on values 2 and 1 in regions 2, 1. 
The result of the integration is 




We have assumed that Vt — c s t > r min} and that the effective size of the background 
medium exceeds Vt + c s t . In the limit M. » 1, we have I sup ersonic — * ln(Vt/r min ); with 
Vt — > r max , this recovers the steady-state result .@ 

In Figure 3, we plot the dynamical friction force as a function of the Mach number for 
several values of c s t/r min . 



3. Discussion 

The main formal result of this paper is the evaluation, in linear perturbation theory, 
of the gravitational drag force on a massive perturber M p moving on a straight-line 
trajectory through an infinite, homogeneous, gaseous medium of density p and sound 
speed c s . Together with the definition in equation fljjD, the expressions ( ft!]) and (|15|) give 



4 As pointed out by the referee Scott Tremaine, the notion that the maximum impact parameter increases 
as Vt was earlier introduced by J. Ostriker & Davidsen (1968) in a time-dependent analysis of stellar 
relaxation. 
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Fig. 3. — Solid curves: Dynamical friction force in a gaseous medium as a function of Mach 
number Ai = V/c s . Curves correspond to \n(c s t / r m i n ) = 4, 6, 8, 16. Dashed curves: 
Corresponding dynamical friction force in a collisionless medium with particle velocity 
dispersion a = c s and r max =Vt — M.c s t. 
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the dynamical friction (df) drag forces on subsonic and supersonic perturbers, respectively. 
Figure 3 presents the same results graphically, showing how the drag force varies with the 
Mach number Ai of the perturber, and the time over which the perturber has been moving 
with fixed speed V = c s hA. 

For comparison, we have also included in Figure 3 the result for the gravitational drag 
on a particle of mass M p moving through a collisionless medium with the same density p 
as the gaseous medium we have considered, and with a Maxwellian distribution of particle 
velocities with a = c s . From equation (7-18) of |Binney fc Tremaine (1987)| , the collisionless 
df drag force is given by expression (O) with 



\ n f r -^L) [erf(X) - ^Le~ x \ (16) 

V T m in J \/7T 



collisionless 

where X = V/(a\/2). From Fi gure 3, it is clear that (a) for Ai » 1, the collisionless and 
gaseous df forces are identical (as has been previously noted); (b) for M. < 1, the drag 
force is generally larger in a collisionless medium than in a gaseous medium, because in the 
latter case pressure forces create a symmetric distribution in the background medium in 
the vicinity of the perturber; (c) the functional form of the gaseous df drag is much more 
sharply peaked near M. = 1 than it is for the collisionless df drag; perturbers moving at 
speeds near Mach 1 resonantly interact with the pressure waves that they launch in the 
background medium; (d) for a given value of ln(A) = ln(r max /r min ), the peak value of the 
gaseous df force is much larger than the corresponding peak value of the collisionless df 
force; at Ai ~ 1, there is a factor of four difference in the force between the two cases. 
Below, we explore some potential consequences of these results in a variety of astronomical 
systems. 

As a consequence of the stronger gaseous df force than collisionless df drag force for 
supersonic motion, massive objects may make their way more rapidly to the center of a 
star cluster or galaxy if they arrive at the outer edge before, rather than after, the gas is 
turned into stars. For a particular example, we consider the decay of a massive perturber's 
near-circular orbit in a spherical density distribution with a singular isothermal sphere 
profile, p{r) = c 2 s /{2iiGr 2 ) so M(r) = 2c 2 s r/G (here c s respectively denotes the sound 
speed or velocity dispersion for a gaseous or stellar distribution). For this density profile, 
the circular speed is constant, V = c s a/2 . By equating the rate of decrease of angular 
momentum d(M p Vr)/dt to the torque = rF d f, one finds that the time for the perturber's 
orbit to decay from r init to r << r init is given approximately by 



td}{r ini t) M{r init ) 



torbinnu) 4nM„ In (-*— ) 



;i7) 
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for gaseous df. This time is a factor of 0.428 shorter than the corresponding decay time 
under stellar df (cf. |Binney fc Tremaine (198"7]| , eq. 7-25). For galactic applications, the 
implication is that condensed objects that form early (e.g. globular clusters) could spiral 
into the galactic center from a factor 1.5 further in a galaxy than would be predicted by 
stellar df theory (|Tremaine, Ostriker, fc Spitzer (19 75 )|) , within the epoch over which the 
baryonic distribution remains gaseous. For star cluster applications, the shorter df time for 
gaseous distributions may help explain why young, embedded stellar clusters like the Orion 
Nebular Cluster ( [Hillenbrand (1997)| , Hillenbrand fc Hartmann ( 1998 )| ) appear significantly 
more relaxed than expected from stellar df alone; e.g., for the ONC, n-body simulations 
show that stellar df would require a time a factor 3-4 longer than the best estimate of 
the cluster age to establish the observed mass segregation ( [Marshall, Ostriker, fc Teuben 

rrgggi. 



As mentioned above (see Fig. 3), the gaseous df force is strongly depressed for subsonic 
perturbers. For modeling the global evolution of combined star-gas systems in which the 
particle velocity dispersion and gas sound speed are comparable, the strong cutoff of F d f 
for subsonic perturbers implies that setting F d f = for V < c s (e.g. [Saiyadpour, Deiss, fc 



Kegel (1997)|) should yield realistic results. However, in other circumstances it is interesting 



to enquire how the small, but nonzero, df drag on subsonic perturbers (computed in this 
paper) can affect their orbit evolution. 

As a model problem, we consider the decay of a subsonically-moving mass M p on a 
near-circular orbit embedded within a constant-density gaseous sphere of radius r$. For 



constant background density, the angular orbit frequency f2 = y4irGp /3 and orbital period 
tori, — 2ir/fl are independent of the distance r from the center. If this const ant- density 
region represents the core of a nonsingular isothermal sphere with core radius ro and sound 
speed c s , then r Q/c s = y/3; thus, the Mach number for a circular orbit at r is M = y/3r/ro. 
Just as above for decay in a power-law density profile, we can compute the time for the 
perturber's orbit to decay from r init = r /\/3 to r/ as 

t{r f ) M core ri MdM 

t orb 7t&l 2 M p JM f =^L I suhsonic (MY 1 j 

the numerical coefficient equals 0.0204. 

Figure 4a shows how the decay of the orbital radius depends on time and on the mass 
of the perturber relative to the whole core. In Figure 4b, we verify that the assumption of 
a near-circular orbit is valid provided M p /M core << 1, since 

Vr_ = 3 5/2 M p I subsonic (M) 

v lD ~ 2M rnrP _ M s ' 1 j 
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Fig. 4. — (a) Decay in time of radial distance r of massive perturber M p on near-circular 
orbit about the center of uniform- density core of mass M core , radius r , and sound speed 
c s = (GM core /(3r )) 1 ^ 2 ; (b) Radial-to-azimuthal velocity ratio for same situation as in (a). 



-14- 




Fig. 5. — For same situation as in Fig. 4, perturber's radial distance from the center at time 
when the forward wave disturbance has propagated ahead of the perturber by r (see text). 
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which has the limit 2.598M p /M core for M. << 1. In arriving at the results shown in Fig. 4, 
we have used equation ([14]). Its validity depends, however, on the size of the uniform- density 
core exceeding that of the perturbed-density region. This assumption must fail, and the 
df drag force consequently decrease, when the forward wave defining the disturbed-density 
region (cf. Fig. 1) reaches a distance ~ r ahead of the perturber.0 The position r(tf w ) 
of the perturber at this time is plotted as a function of M p /M core in Figure 5. Based on 
this figure, only perturbers of mass M p > 0.2M core are predicted to reach within one-tenth 
of a core radius before the df drag decreases. For lower-mass perturbers, this implies that 
the decay of orbits to very small radii may stall, if the df drag is sharply reduced after the 
time tf w . This result may have relevance for models of QSO evolution (|Silk fc Rees (1998)| ) 



in which primordial black holes are formed away from the centers of galaxies, later to be 
driven there by df during mergers. If, as assumed in this scenario, these events occur before 
the advent of star formation, than the relevant df drag is the gaseous df examined in this 
paper. If the df drag becomes inefficient when the orbit reaches the protogalaxy core and 
becomes subsonic, then these massive black holes may have more time to grow by accretion 
before they finally sink to the centers of their host galaxies. 

Finally, we comment on the applicability of our results to the interaction of 
protoplanets with a gaseous circumstellar nebula in which they may grow. Protoplanets 
orbit faster than the surrounding gaseous nebula, with the Mach number of the drift speed 
M. = v re i/c s ~ c s /(flK r ) « 1 (where VL K is the local Kepler angular rotation frequency 
and r is the protoplanet's semimajor axis). Drag forces occur as a result of this relative 
motion, leading to inward radial migration of protoplanets as they lose angular momentum 
to the surrounding nebula (e.g. [Ward (1997)| and references therein). Using parameters 
based on solid bodies in the protosolar nebula, it can be shown that the nominal df drag 
from equations (R) and (0), F df « (Air/3)(GM p ) 2 p /(c s n K r) « (2tt /3)(GM p ) 2 i:/(c 2 r), 



exceeds the large- Reynolds number aerodynamic drag F aer « {l/4)v 2 el iiR 2 pQ ( |Landau ~EZ 
Lifshitz (1987)| , §45) for protoplanets of radius R p greater than a few 100km. 



It is sometime argued that gravitationally-enhanced drag on protoplanets (cf. Takeda 



et al (1985)| , phtsuki, Nakagawa, fc Nakazawa ( 1988 )| , but note the different scaling from our 



results) may enhance migration rates over those predicted to arise from differential resonant 



5 Because the perturber follows a circular rather than straight-line orbit, the df drag at late times should 
still be nonzero: Since the direction of V changes by tt/2 four times per orbit, the forward- wave propagation 
effectively "restarts" as well. When r << ro, over each quarter orbit, each new expansion wave propagates 
to just ~ r ; over much of the quarter-orbit, an unbalanced trailing density enhancement will remain within 
the core. Thus, we expect that the df drag will still be some fraction of the value found using eq. (|l4|), 
although a more refined calculation is needed to predict the fraction quantitatively. 
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torques flGoldrcich fc Tremaine ( 1980 )|) between the protoplanet and the surrounding 
nebula. In fact, provided that the difference of gas orbits from Kepler orbits is included so 
that the perturber does not corotate with the nebular gas at the same radial distance, then 
the df drag is automatically incorporated in the net resonant torque. Indeed, calculations 
show that the net resonant torque from a thin disk (e.g. [Korycansky fc Pollack (1993)1 , [Ward 



( 1997)|) , taking this velocity difference into account, is of the same magnitude (and scaling) 
as the df drag estimated above. Thus, while the gaseous df drag may have implications 
for the late stages of planet formation, its effects are naturally incorporated within a full 
resonant formalism. Although such a calculation has not yet been performed in three 
dimensions, it is reassuring that our simple estimate of the drag - which neglects gradients 
of velocity, density, and temperature in the disk - and the estimates from 2D resonant 
torque differences - which neglect the disk thickness - nevertheless yield comparable 
answers. 

The author is grateful to the referee Scott Tremaine for helpful comments on the 
manuscript. 
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